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Abstract 

This work presents theoretical calculations of the electrical transport properties of the Ag, Au, 
and La fractionally filled bulk skutterudites: CoSb3, C0AS3, and C0P3. Density functional theory 
(DFT), along with projector augmented wave (PAW) potentials, was used to calculate bulk band 
structures and partial density of states. The Seebeck coefficient (S), electrical conductivity (a), 
and power factor (S 2 a) were calculated as a function of temperature and filling fraction using 
the momentum matrix method along the entire first Brillouin zone. Calculated trends in the 
electrical transport properties agree with previously published experimental measurements for p- 
type unfilled and La filled CoSb3. The calculated S and a values for the Ag and Au filled compounds 
indicates that the most promising electronic properties are exhibited by p-type Auo.i2s(CoSb3)4, 
Auo.25(CoSb3)4, and Au(CoSb3)4. Au is therefore recommended as a promising filler for improved 
thermoelectric properties in cobalt antimonides. Ag is also a good filler for cobalt phosphides; the 
creation of a negative indirect band gap is observed in Ag(CoPs)4, which indicates semimetallic 
behavior, so this compound may possibly exhibit lower thermal conductivity than metallic C0P3. 
Finally, we recommend future directions for improving the thermoelectric figure of merit of these 
materials. 

PACS numbers: 72.20.Pa, 71.20.Nr, 71.15.Mb 

Keywords: density functional theory, filled skutterudite, thermoelectric power factor, electrical transport 
properties 



I. INTRODUCTION 



Thermoelectric materials can be used for solid state power generation in a variety of appli- 
cations, including automotive power generation 1 ' 2 and solar energy capture and conversion 3 . 
Such devices operate without the need for moving parts, and can be incorporated into most 
current energy generation technologies to capture waste heat and increase conversion effi- 
ciency Current thermoelectric devices are relatively inefficient (~ 18% 4 ) but with recent 
advances in nanostructuring concepts, it has been shown that more efficient devices may be 
attainable 5 . 

In this study, we focus specifically on skutterudite materials, which are attractive for 
thermoelectric applications because of their large unit cells and high intrinsic electrical 
properties. The skutterudite unit cell consists of eight TM 6 octahedral groups - one in 
each corner of the cubic unit cell - where T is a transition metal (e.g., Co, Rh, or Ir) and 
M is a metalloid or nonmetal (e.g., Sb, As, or P) (Figure 1). Each M atom belongs to 
two adjacent octahedra and each T atom is at the center of the octahedral group. These 
octahedra are arranged in such a way that a cage-like void is created in the center of the 
unit cell. These voids are large enough to incorporate certain atoms (X) that, when present, 
have the ability to lower the thermal conductivity of the material by "rattling" about within 
the void, decreasing the phonon mean free path, and thus, increasing phonon dispersion. 

Previous studies have shown that filling the naturally occurring voids within the skut- 
terudite lattice with heavy element atoms, (e.g., Lanthanides La, Ce, and Eu), leads to 
a considerable reduction in the thermal conductivity 6 ' 7 . In general, ideal filler atoms are 
described as small and heavy, with electrons that can readily hybridize with the conducting 
electrons in the lattice 8,9 . These hybridized electrons may help prevent large reductions in 
electronic properties that might otherwise occur through void filling. Recent studies have 
indicated that the filling of the lattice voids with smaller or lighter atoms with more metallic 
properties, such as Ba 10 , Na 11 or In 12 , leads to an increase in the thermoelectric figure of 
merit compared to that obtained by filling with La. This is in part attributable to the en- 
hancement of electronic properties from metallic fillers as compared to heavier Lanthanides. 
Furthermore, it has been shown that double element filling may further enhance thermo- 
electric performance 13 compared to single element filling, though this is outside the scope of 
the present work. 
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In this study, we test the accuracy of density functional theory (DFT) calculations, fol- 
lowing the Projector Augmented Wave (PAW) approach, for predicting electrical transport 
properties for skutterudite thermoelectric materials; we use the momentum matrix method 
in conjunction with PAW (PAW-MM) and compare our calculated results to previously pub- 
lished experimental results. Many previous theoretical investigations into skutterudites have 
utilized the full-potential linearized augmented plane wave (FLAPW) approach to calculate 
electronic properties, but we are specifically interested in PAW for its computational speed 
and efficiency. These qualities are desirable when screening a large group of prospective 
materials for a specific application. We calculate the thermoelectric properties for both 
well-characterized and previously uncharacterized skutterudites. We emphasize that the 
computational method employed is a simplified calculation for ideal bulk crystal systems 
that does not take into account grain boundary effects or the dependence of the electronic 
time relaxation constant, r e , on temperature. Therefore, only trends and relative magnitudes 
of properties are demonstrated here. 

We also aim to determine the effect of filler size, atomic weight, and electronegativity on 
the resulting electrical transport properties of filled Co-based skutterudites. A summary of 
filler properties is presented in Table I. We hypothesize that Au is a promising metallic filler 
atom due to its small size and large mass, and because its valence electron configuration is 
similar to that of La, which has been previously shown to be a successful filler. Both Au and 
La have similar core electronic structures and valence shells composed of all but one filled 
orbitals and one unpaired electron in the highest energy orbital; they differ only in that the 
unpaired electron in La is a d electron and that in Au is an s electron. Furthermore, Au is 
heavier (atomic weight is 197 amu) and smaller (atomic radius is 135 pm) than La (atomic 
weight is 139 amu, atomic radius is 195 pm). Both of these desirable characteristics should 
result in improved thermoelectric properties for Au filled skutterudites. We also study Ag as 
a candidate filler, since Ag has a similar electronic structure compared to Au, but is lighter 
(atomic weight is 108 amu) and larger (atomic radius is 144 pm). From size considerations, 
we would expect that the size of the filler atom to be important if the atom is large enough 
to interact with the void boundaries; otherwise, the heaviness of the atom is likely to be 
more important for phonon dispersion. Finally, we note that the electronegativities of Ag, 
Au, and La are 1.9, 2.4, and 1.1, respectively, and those of P, As, and Sb are 2.1, 2.0, and 
1.9, respectively. We expect that larger differences in electronegativity between the filler 
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and host (M) atoms would improve electron transport in filled skutterudites. All of these 
effects will be explored in the calculation of the thermoelectric strength of Ag, Au, and La 
filled Co-based skutterudites. 

II. METHODS 
A. Theory 

One quantity that is commonly used to characterize the thermoelectric strength of a 
material is the thermoelectric figure of merit 4 , z: 

zT= T=— . 1 

k pn 

Here, S is the Seebeck coefficient, a is the electrical conductivity (with p being its inverse, 
or the electrical resistivity), and k is the thermal conductivity. We aim to maximize S and a 
and minimize k for optimal thermoelectric performance. Equation 1 suggests that for every 
increase in the Seebeck coefficient, AS, there is a corresponding increase, (AS), in zT. a 
and k are directly and inversely proportional to zT, respectively. In this work, we will focus 
on optimizing the electrical transport properties of skutterudites by characterizing S and a. 

Traditionally, band structure calculations have been used to directly calculate electronic 
properties. However, the problem with this approach is that many DFT software packages 
output the band structure in the form of an array organized by increasing energy, without 
band distinction. Thus, errors may arise in the computation of the electrical conductivity, 
particularly at band crossings around the Fermi level, even when Fourier transformations 
are used 14 . A different approach is desirable in complex systems such as skutterudites, when 
numerous band crossings are present. 

The momentum matrix method rectifies this issue by relying directly on the wavefunction 
descriptors. We can express the computed electron group velocities as follows 15 : 

v a f} (i, k) = — (i, k) \p\ ij af} (i, k)) (2) 

Here, v are the velocities, ip are the wavefunctions, m is the electron mass, and p is the 
momentum operator, {i, k) denotes the eigenvalue-band index and a(3 denote the directional 
coordinates. The resulting matrix of velocities indexed by energy and /c-point is called the 
optical matrix. Although this equation may be used to determine group velocities from exact 
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wavefunction descriptions, the PAW-MM approach requires the addition of an augmentation 
factor to account for errors arising from the pseudization of the wavefunctions around the 
atomic cores 16 ' 17 . 

Boltzmann transport theory is then used to derive the following nested equations that 
enable us to determine a and S in terms of the group velocities 14 : 

(i, k) = e 2 r i:k v a (i, k) vp (i, k) (3) 

^(0 = ^E^(a) * (e ~ e ^ (4) 

i,k 
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Here, N is the number of Appoints, 5 is the unit impulse function, Q is the volume of the 
primitive cell, ep is the Fermi level, and / is the Fermi-Dirac distribution for an electron 

gas: 
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Equations 5-7 yield second-rank tensors. To quantify these properties as scalars, the trace 
of each matrix is evaluated. 



B. Calculations 



The electrical transport properties of fractionally filled Co-based skutterudites - specifi- 
cally, CoSb 3 , CoAs 3 , and CoP 3 with Ag, Au, and La filler atoms - were evaluated. La was 
chosen for its previously demonstrated ability to increase the thermoelectric figure of merit 
of filled skutterudites; thus, La will serve as a baseline for comparison to the previously 
uncharacterized Ag and Au fillers. 

All calculations of electronic wavefunctions were performed using DFT, as implemented 
in the Vienna Ab Initio Simulation Package (VASP) 18-21 . PAW potentials 22 ' 23 were used to 
maintain a balance between high accuracy and reasonable computational costs. The general- 
ized gradient approximation of Perdew, Burke, and Ernzerhof 24,25 was used to represent the 
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exchange-correlation energy. All calculations were performed using the blocked Davidson it- 
erative matrix diagonalization algorithm 26 to optimize the orbitals. The tetrahedron method 
with Blochl corrections was used to determine partial occupancies. The first Brillouin zone 
was sampled with 100 fc-points, and the kinetic energy cutoff for the plane wave basis set 
was set to 270eV after extensive fc-point convergence tests. The resulting wavefunctions are 
used in the momentum matrix equation 2, with a PAW potential correction to determine 
the group velocities as noted in Section II A. Equations 3-7 are then used to determine the 
thermoelectric properties from the group velocities with a code developed in MATLAB 27 . 

As shown in Figure 1 there are two voids per unit cell, so that only a limited number 
of filling fractions (x G {0.5, 1}) may be explored. In order to obtain a larger sampling 
set of filling fractions, a 2 x 2 x 2 supercell containing 16 voids was constructed. Voids 
were filled in pairs until the structure was fully filled. Table II shows the resulting highest 
symmetries used for the filled systems, X x (CoM 3 ) 4 . S, a, and S 2 a values were calculated 
for x G {0.125,0.25,0.375,0.5,0.625,0.75,0.875,1}; a denser sampling of filling fractions 
would be difficult to achieve given computational limitations and may not yield significantly 
different trends in the resulting data. All calculations were performed at the following 
temperatures: T G {300, 400, 500, 600, 700, 800, 900} K. 

Calculations were first performed on the unfilled CoSb 3 system and the La filled, Fe doped 
CoSb 3 system at different carrier concentrations, and compared to previously published 
experimental data. Doping was mimicked using the rigid band model (i.e., moving the Fermi 
energy on a fixed band structure) to simulate different electron carrier concentrations, n e , 
and hole carrier concentrations, n^. These carrier concentrations were determined using the 
density of states, g (e): 



1 oo 

n e = v Y,9c (e) f (e) As (9) 

* e=ep 

nu= \ £ 9v(e)(l-f(e))Ae (10) 

The c and v subscripts denote the conduction and valence states respectively. Then, the 
Ag, Au, and La filled (undoped) skutterudites were screened across the range of x and T 
values in the sample space, and the electrical transport properties - S, a, and S 2 a - were 
calculated from the wave descriptors. Additional calculations of the electrical transport 
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properties were performed on the most promising filled systems for a range of deviations in 
the chemical potential. 

III. RESULTS AND DISCUSSION 

A. Unfilled skutterudites 

The rightmost column in Table III shows the calculated void radii as determined upon 
geometry optimization of the skutterudite systems. Full lattice optimizations resulted in 
only a +0.6% increase in the lattice parameters; thus, we allowed only the atoms to relax 
during geometry optimization. Smaller voids are observed in C0AS3 and CoSb 3 and larger 
voids are observed in C0P3 upon geometry optimization. 

Figure 2 shows the band structures calculated for the three unfilled skutterudites (i.e., 
C0P3, C0AS3, and CoSb 3 ). CoSb3 is a semiconductor with a small band gap of 0.3 eV, 
C0AS3 has a pseudo gap of eV, and C0P3 is a conductor with an overlap at the T point; 
these results generally agree with previous findings 28 ' 29 . The total density of states (DOS), 
presented alongside the partial contributions from the Co 3d electrons and the 3p, 4p, and 
5p electrons of P, As, and Sb, respectively, are also shown in Figure 2. The dotted horizontal 
lines denote the Fermi-Dirac distribution filter (§f ) of equations 5-6 at 300K in this and all 
subsequent figures. 

From these plots we observe that there is significant hybridization of the valence orbitals 
at the Fermi level for all three bulk materials. We also notice that a majority of the DOS 
contribution around the Fermi energy may be attributed to the valence electrons. These DOS 
and band structure diagrams indicate that the valence electrons are essential in dictating the 
electronic properties of the compound, and fillers containing electrons that may hybridize 
with electrons at these energy levels can strongly impact electrical transport properties. 

B. Comparison of Computational to Experimental Results 

We first compare the calculated Seebeck coefficient, S, and electrical resistivity, p, to 
available experimental values for two systems - unfilled CoSb 3 and La filled, Fe doped 
CoSb 3 . Figure 3 shows that the trends in both properties and the peak in S are predicted 
correctly for unfilled p-type CoSb 3 at high hole carrier concentrations. By comparing our 
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calculated values of p to the experimental data at T = 500 K, we determine that the best-fit 
value for the electronic time relaxation constant for unfilled CoSb 3 is r e = 5 x 10~ 13 s. For 
low carrier concentrations, (i.e., where the Fermi energy approaches the undoped value), 
our theoretical values of S drop off more rapidly than the observed experimental values. At 
room temperature, our calculated S values are correct; however, as the temperature rises, 
the discrepancy increases. Since we determine our band occupations from the Fermi dis- 
tribution, it is possible that the deviations arise from the inclusion of conduction bands in 
the calculation. As the temperature increases, more conduction bands are included in the 
computation. Furthermore, as the hole carrier concentration increases {ep moves downward 
on the band structure), the Fermi distribution includes fewer conduction electrons and our 
estimation improves. This is a reflection of the limitations of DFT in describing conduc- 
tion band characteristics, especially at band edges. Although the theoretically predicted 
and experimentally measured values of p and S do not always agree throughout the sam- 
pled temperature range, the trends in these properties with temperature are consistent and 
demonstrate that there is an optimal temperature range for high values of S and low values 
of p. 

Although not directly the focus on this work, we note that that there have been prior 
comparisons made between theory and experiment for unfilled n-type CoSb 3 30,31 , so we now 
demonstrate that our model is flexible enough to consider both types of dopants. For n-type 
materials, the number of bands included in the calculation must be increased in order to 
ensure convergence of the behavior of the conduction bands. Figure 4 shows a comparison 
between the experimental trends and the theoretical trends for unfilled n-type CoSb3. For 
high electron carrier concentrations, where the Fermi energy lies within the conduction band, 
we achieve good agreement between theory and experiment. The carrier concentration used 
in the theoretical calculations was obtained by fitting the computed S to the experimental 
data at T = 300 K. As depicted in Figure 4, the shape of the computed curves does not 
exactly match the experimental data across the range of temperatures sampled, although the 
computed values of S at low temperatures match within 10% of the experimental values. We 
attribute this discrepancy to the observation that in experimentally developed materials, the 
carrier concentration and transport properties will fluctuate differently with temperature, 
as a result of defects, than what would be expected in a defect-free model like that used in 
our calculations. Deviations in the steepness of the curvature of the bands are thus most 
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likely due to the constant relaxation time estimation 31 used in our calculations. For very low 
carrier concentrations (not shown), where the chemical potential falls within the band gap, 
we observe a greater deviation from the experimental results, especially at low temperatures. 
Since the Fermi-Dirac distribution does not accurately describe the distribution of electrons 
within the band gap, we must instead increase the k-point sampling of the Brillouin zone to 
5,000 points and triple the number of total bands in the calculation to achieve reasonable 
accuracy within the band gap. However, we again stress that for the filled systems considered 
subsequently in this study, the filler acts as a p-type dopant and the Fermi level is pushed 
into the valence band, so that only heavily doped materials are observed; therefore, such 
resource-intensive calculations on n-type systems are outside the scope of the present work. 

The results for La filled, Fe doped skutterudites are shown in Figure 5. The experimental 
measurements are made on 90% La filled skutterudites, but given the range of filling fractions 
that were computationally tractable, we relate these results to our calculated predictions on 
100% La filled skutterudites. We see from these plots that the trace of the Seebeck coefficient 
has a flatter temperature dependence compared to the experimental trends. However, after 
investigating the directional components of S, we see that the S yy component exhibits a 
similar trend to the experimental values. As we adjust the value of n to attempt to emulate 
90% La filling, we notice that the values asymptotically approach the experimental results. 
Thus we conclude that the small discrepancy between the calculated and experimental values 
of S is likely due to the difference in the La filling fraction. Again, we fit the calculated 
values of p to the experimental data at T = 500 K, and determine that the best-fit value 
for the electronic time relaxation constant for La filled, Fe doped CoSb 3 is r e = 2 x 10" 13 s. 
From these results, it seems reasonable to estimate the electronic relaxation time constant 
for all of the systems considered in this study to be r e = 10~ 13 s, and we will use this value 
for all subsequent calculations. 

C. Ag and Au filled skutterudites 

We next explored Ag and Au filled skutterudites as potential candidates for improved 
thermoelectric behavior. We plotted the calculated power factor, S 2 <r, which is the numerator 
of the figure of merit, as a function of temperature. For these calculations, the time relaxation 
constant was again assumed to be independent of filling fraction and temperature, and was 
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approximated as r e = 10~ 13 s, as described in Section III B. Figure 6 shows the power factors 
plotted as a function of temperature for the undoped compounds at the filling fractions that 
yielded Im3 symmetry of the unit cell. Charge compensation was not included so that the 
effect of composition on the thermoelectric power factor may be isolated and identified as 
follows: 

• C0P3: Ag is a reasonably good filler at 100% filling fraction above room temperature. 
Au, however, is not a good filler in this compound. La is an excellent filler and actually 
enhances electronic properties for fully filled lattices. 

• C0AS3: Both Au and Ag are good fillers at temperatures above 500 K. The power 
factor for Ag(CoAs 3 ) 4 experiences a dip at 450 K because of the transition from p- 
type to n-type behavior, as reflected by the sign of S. La is a moderately good filler. 
However, completely filling the lattice with Au seems to have an enhanced effect on 
electrical transport properties compared to partial filling, so this system is worthwhile 
as the subject of future research for high temperature applications. 

• CoSb 3 : Ag is a poor filler. Au holds great promise as a filler, especially for low 
filling fractions at temperatures between 500 K and 800 K and for complete filling 
at temperatures above 700 K. Indeed, complete filling with Au exhibits the greatest 
enhancement of electrical transport properties at very high temperatures, even when 
compared with La, which is an excellent filler in fully filled systems. 

Figure 6 indicates the potential of fillers for high-performing thermoelectric behavior, 
based on undoped systems. The power factors at 100% filling are markedly higher because of 
a multiple- fold increase in the electrical conductivity, specifically in the a zz component. This 
is likely due to increased coordination between the filler atoms and the host lattice in the fully 
filled configuration. Ag and Au are excellent conductors, and in a fully filled skutterudite 
lattice, the layout of the fillers is similar to that in the pure conductors. Therefore, when the 
voids are fully filled, there may be resonance at the location of the fillers, so that a highly 
electrically conductive material would be produced. Experimentally, this is corroborated by 
90% La filled CoSb 3 exhibiting the best performance among similar skutterudites 6 . 

Although the fully filled systems in general show the greatest promise, 100% filling is 
rarely achievable without heavy doping for charge compensation. Furthermore, a system 
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that does not possess high power factors when undoped may still be a good thermoelectric 
when doped, since it is also likely that fully filled materials would exhibit high thermal 
conductivities without the introduction of impurities or defects. In any event, screening cal- 
culations on undoped systems, as performed here, can aid in highlighting systems where the 
filler acts as a dopant and enhances electrical transport properties. For the examples of solar 
thermal or powering electronics at mid-range temperatures (e.g., 500-600 K), fractionally 
filled materials, such as the Au .i25(CoSb 3 ) 4 system, show the most promise for improved 
thermoelectric behavior. 

We note that Ag performs well only for the C0P3 system, while Au performs well only 
for the C0AS3 and CoSb3 systems. Thus, it appears that filler size is important only when it 
exceeds a threshold percentage of the void radius. Since Au is the smallest filler considered 
in this study, it does not fulfill this criterion for any Co-based skutterudite. La, being the 
largest filler considered in this study, does fulfill this criterion for all Co-based skutterudites 
and exhibits good performance in all cases. However, higher filler weights are important 
when considering systems with larger voids, and will likely have an even greater impact on 
thermal behavior. 

The electronegativity of the filler may play an even greater role in determining the elec- 
trical behavior of the system. La, with an electronegativity of 1.1 on the Pauling scale, 
readily donates electrons to the lattice, which results in enhanced electrical properties. On 
the other hand, Au and Ag atoms are more selective in their willingness to donate electrons. 
Au has an electronegativity of 2.4 while Ag has an electronegativity of 1.9 on the Pauling 
scale. The electronegativity of Ag is much closer to that of the M atoms in the skutterudite 
(1.9-2.1 on the Pauling scale), while the difference in electronegativity between Au and the 
M atoms is much larger. In fact, since the electronegativity of Ag is so close to that of the 
M atoms, it is a poor filler in all but the phosphide skutterudites. On the other hand, Au is 
a good filler in arsenide and antimonide skutterudites, where the electronegativity difference 
is 0.4 and 0.5, respectively, but poor in the phosphides. Thus, we propose that the absolute 
difference in electronegativity between the filler and host atoms, combined with the atomic 
weight of the filler atom, are important considerations for designing improved thermoelectric 
materials. 

We now examine in greater depth four systems that exhibit improved electrical transport 
properties upon filling: 
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• 12.5% Au filled CoSb 3 

• 100% Au filled CoSb 3 

• 100% Au filled CoAs 3 

• 100% Ag filled CoP 3 

Figure 7 shows the band structures and density of states for these systems. Figures 8-9 
show the computed values of the Seebeck coefficient, S, and electrical conductivity, a, as a 
function of temperature for these systems. Two observations may be made: 

• From the density of states, it seems that there is no appreciable band gap in the fully 
filled skutterudites, Au(CoSb 3 ) 4 , Au(CoAs 3 ) 4 , and Ag(CoP 3 ) 4 , so it would appear 
likely that these compounds would be conducting and would not make good ther- 
moelectric materials. However, in the case of Au(CoAs 3 ) 4 and Ag(CoP 3 ) 4 , the band 
structure indicates that the filler creates an indirect band gap of -0.8 eV and -0.45 
eV, respectively; these indirect band gaps are seen between the P and T Appoints. We 
expect that these materials would therefore exhibit semimetallic rather than metallic 
properties. The sign of the S values indicate that Au(CoAs 3 ) 4 exhibits p-type behav- 
ior while Ag(CoP 3 ) 4 exhibits n-type behavior. These materials should thus be further 
studied to ascertain their thermal conductivities. 

• Filling CoSb 3 with small fractions of Au does not appreciably alter the narrow band 
gap. In Au .i25(CoSb 3 ) 4 , S is greatly enhanced while a experiences only a slight 
increase. The resulting undoped, Au filled compound is a p-type semiconductor with 
high S values. To improve the values of a, further analysis of doping was performed 
to determine whether enhancement of electrical transport properties is possible even 
at low Au filling fractions; this will be described in Section III C. 

Figure 10 shows a complete thermodynamics cycle for the formation of the minimally 
filled Au .i25CoSb 3 system from CoSb 3 . 

• First, we calculated the total energy of unfilled CoSb 3 (-646.6 eV), and set this to be 
our reference point. 

• Second, we calculated the total energy of 12.5% Au filled CoSb 3 (-648.9 eV). 
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• Third, we calculated the total energy of unfilled CoSb 3 with the same void radius 
as 12.5% Au filled CoSb 3 , by simply removing the Au filler atom and performing an 
energy minimization without geometry optimization. 

By subtracting the total energy of the Au filler atom (-0.1 eV) from the system, we see that 
the formation of the filled compound is thermodynamically favorable, with an overall energy 
of formation of -2.2 eV. Thus, it should be theoretically possible to synthesize partially Au 
filled skutterudites and experimentally characterize its properties for use in thermoelectric 
devices. 

We also performed an analysis of the maximum filling fraction possible based on previ- 
ously established methods 32 . This analysis tells us that, depending on the structure of the 
alternate stable phase AuSb 2 , a maximum filling fraction of anywhere between 6% (Fm3m 
symmetry) to 100% (P2iC symmetry) is possible without doping. Further information about 
secondary phases is necessary to more accurately determine the maximum filling fraction. 

D. Analysis of unit cell symmetry 

We previously presented our calculated electrical transport properties for the filled skut- 
terudite systems exhibiting Im3 symmetry in their unit cells, so we had omitted the 25%, 
50%, and 75% filling fractions in Figure 6. For Au^CoSba)^ the 25% filling fraction has 
three possible symmetry configurations: R3, Cmmm, and Pmmm, the 50% filling fraction 
has five possible symmetries, and the 75% filling fraction has three possible symmetries. 
Since the 12.5% filling fraction had been shown in Section III C to exhibit improved elec- 
trical transport properties, we were especially interested in exploring all of the possible 
symmetries for the 25% filled structure, but in the case of the 50% and 75% filled struc- 
tures, we chose instead to perform calculations for only one symmetry group each, which 
were Pm3 and Cmmm, respectively. We note that Im3 is considered a high symmetry group, 
so the other symmetry groups being considered are lower in symmetry than Im3 (Figure 11). 
In Figure 12, we plotted the power factor, S 2 a, of CoSb 3 as a function of Au filling fraction 
at different unit cell symmetries corresponding to the filling fraction. We observe that visible 
distortions in the overall trends are indeed observed at the 25%, 50%, and 75% Au filling 
fractions, compared to those computed with Im3 symmetry. 

To quantify the effect of symmetry on electrical transport properties, the 25% Au filled 
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system was re-optimized at all of its possible symmetries - Cmmm (orthorhombic dipyra- 
midal), Pmmm (orthorhombic dipyramidal), and R3 (trigonal rhombohedral). From Figure 
12, we see that the R3 symmetry for the unit cell results in the largest power factors. We 
conclude that the pattern in which the voids are filled is important for achieving favorable 
electrical transport properties, but stress that the full range of predicted values is likely to 
be observed experimentally when symmetry is difficult to control. Figure 12 also suggests 
that 25% Au filling at Cmmm symmetry may exhibit a higher power factor than 12.5% Au 
filling at Im3 symmetry, and this prediction is confirmed by the results shown in Figure 13 
for power factor as a function of temperature. 

Figure 14 shows the S, a, and S 2 a values for the 12.5%, 25%, and 100% Au filled CoSb 3 
systems at 500 K. The graphs show values for doping levels from -0.5 eV to +0.5 eV around 
the Fermi energy, as simulated using the rigid band model. From these graphs we note 
that filling with Au creates a second hump in the S value which shifts further away from 
the fermi energy with increasing filler concentration. This second hump also results in the 
highest values of S 2 a. The greatest potential for thermoelectric behavior, in terms of the 
tr (S) value, is observed in the 25% filled compound, while directionally, it is observed for 
the S xx component in the 100% filled compound. Fully filling the lattice with Au appears 
to unfairly dope the host material by a suboptimal amount (i.e., too little or too much) 
into the valence band, so that the power factor is compromised at 500 K. However, further 
improvement in the power factor may be achieved by doping the material. We conclude that 
low fractional filling with Au appears to induce promising p-type behavior and provides a 
good complement to n-type Lanthanum filled materials in p-n junction devices. 



IV. CONCLUSIONS 



We have shown that density functional theory calculations following the Projector- 
Augmented Wave approach, in conjunction with the momentum matrix method for calcu- 
lating electrical transport properties, may be used successfully to obtain predictive trends 
in the performance of Co-based skutterudites over a temperature range of 300-900 K. Au 
filling results in p-type behavior whereas Ag filling results in n-type behavior. When con- 
sidering the properties of candidate fillers, the relative weight and electronegativity appear 
to affect electrical transport properties more than the absolute size of the filler atom. These 
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properties are evident when considering that the hybridization of the filler atom's s elec- 
trons is most emphasized in the fully filled skutterudites, as demonstrated in the density 
of states around the Fermi level in Figure 7. Therefore, in these fully filled systems, a 
resonance quality arises in the lattice where the coordination of electrons results in greatly 
enhanced electrical conductivities; this enhancement is most notable at high temperatures. 
Based on our calculations, fully Au filled CoSb3 would make an excellent thermoelectric 
material at temperatures above 600 K, but it may be difficult to synthesize this material 
in a cost-effective manner 33 ' 34 and it may even exhibit undesirable (semi) metallic behavior, 
as suggested by our band structure calculations. Thus, we propose that low Au filling 
fractions (e.g., 25% or 12.5%) offer the optimal balance between improved performance 
and economical value. In terms of building an efficient thermoelectric device with a p-n 
junction, we suggest that coupling p-type Au filled CoSb 3 with n-type Lanthanum (or other 
rare earth element) filled CoSb 3 would favorably result in a device with optimal performance 
yet minimal lattice strain at the heterojunction. 
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TABLE I. Filler characteristics 



Filler Size (pm) Atomic weight (amu) Electronegativity 
Ag 144 108 1.9 
Au 135 197 2.4 
La 195 139 1.1 



TABLE II. Unfilled (CoM 3 ) and filled (with X) Co-based skutterudites 



Filling Fraction 



Chemical Formula 



Highest Computed Symmetry Space Group 





0.125 
0.25 
0.375 

0.5 
0.625 
0.75 
0.875 
1 



CoM 3 

X .125(CoM 3 )4 

Xo.25(CoM 3 ) 4 

X .375(CoM 3 )4 

Xo. 5 (CoM 3 )4 

X .625(CoM 3 ) 4 
X .75(CoM 3 ) 4 
X . 8 75(CoM 3 ) 4 

X(CoM 3 ) 4 



Im3 
Im3 
R3 
Im3 
Pm3 
Im3 
Cmmm 
Im3 
Im3 



TABLE III. Structural parameters for unfilled skutterudites. 



System 


Lattice Parameter (A) 
Experimental 21 


Void Radius (A) 
Experimental 21 Computed 


CoSb 3 


9.0385 


1.892 


1.875 


CoAs 3 


8.2055 


1.825 


1.780 


C0P3 


7.7073 


1.763 


1.768 



a Experimental values from Nolas et al. 
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FIG. 1. Skutterudite (C0M3) unit cell structure. The top unit cell shows the sixfold coordinated 
Co atoms (blue). The bottom unit cell shows a filler atom (gold) inserted in the center void. The 
second void in the unit cell is located at the corners. 
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FIG. 2. Band structure and partial density of states for unfilled skutterudites. 
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FIG. 3. Comparison of theoretically computed and experimental 36 electrical transport properties, 
S and p, for unfilled p-type CoSt>3, at different hole carrier concentrations, n. 
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FIG. 4. Comparison of theoretically computed and experimental 31 electrical transport property, 
S, for unfilled n-type C0SD3, at different hole carrier concentrations, n. 
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FIG. 5. Comparison of theoretically computed and experimental 6 electrical transport properties, 
S and p, for La(Coo.25Feo. 75803)4 and Lao.g(Coo.25Feo.75Sb3)4, respectively. 
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FIG. 6. Comparison of Ag, Au, and La fillers for C0M3 (M = P, As, Sb) skutterudites at filling 
fractions corresponding to Im3 symmetry in the unit cell. 
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FIG. 7. Band structures and partial density of states for best-performing Ag and Au filled Co- 
based skutterudites: 12.5% and 100% Au filled CoSb 3 , 100% Au filled CoAs 3 , and 100% Ag filled 
CoP 3 . 
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FIG. 8. Electrical transport properties, S, a, and S 2 a, as a function of temperature for 12.5% and 
100% Au filled CoSb 3 and 100% Au filled CoAs 3 . 
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FIG. 9. Electrical transport properties, S, a, and S 2 a, as a function of temperature for 100% 
filled C0P3. 
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FIG. 10. Thermodynamic cycle indicating the thermodynamically downhill energy of formation 
for Au filled CoSb 3 . 
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FIG. 11. Symmetries of filled skuttemdite systems. Cubic systems have the highest symmetry, 
followed by rhombohedral and then orthorhombic systems. 
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FIG. 12. Electrical transport properties, S, a, S 2 a, as a function of Au filling fraction, x, for all 
symmetries computed for Au^CoSbs)^ 
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FIG. 13. Thermoelectric power factor, S 2 a, as a function of temperature for 25% Au filled C0SD3. 
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FIG. 14. Directional (i.e., xx, yy, and zz) electrical transport properties, S, cr, and S ,2 cr, as a 
function of energy for 12.5% (Im3), 25% (R3), and 100% (Im3) Au filled CoSb 3 . 
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